Integration by parts method